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ABSTRACT 


Helioseismic and magnetohydrodynamic waves are abundant in and above 
sunspots. Through cross-correlating oscillation signals in the photosphere observed by 
the SDO/HMI, we reconstruct how waves propagate away from virtual wave sources 
located inside a sunspot. In addition to the usual helioseismic wave, a fast-moving 
wave is detected traveling along the sunspot’s radial direction from the umbra to 
about 15 Mm beyond the sunspot boundary. The wave has a frequency range of 
2.5 —4.0 mHz with a phase velocity of 45.3 km s _1 , substantially faster than the typ¬ 
ical speeds of Alfven and magnetoacoustic waves in the photosphere. The observed 
phenomenon is consistent with a scenario of that a magnetoacoustic wave is excited 
at approximately 5 Mm beneath the sunspot, and its wavefront travels to and sweeps 
across the photosphere with a speed higher than the local magnetoacoustic speed. The 
fast-moving wave, if truly excited beneath the sunspot’s surface, will help open a new 
window to study the internal structure and dynamics of sunspots. 

Subject headings: Sun: helioseismology — sunspots — Sun: oscillations 


1. Introduction 


Oscillations, helioseismic waves, and different types of magnetohydrodynamic (MHD) waves 
are abundant in sunspots’ atmosphere, both at and above the photospheric level (Bogdan & Judge 
20061) . Studying these various waves can help us understand the physical conditions, such as tem¬ 
perature, density stratification, and magnetic field structure, in the atmosphere and even beneath the 
photosphere of sunspots. Inside sunspots, 5-minute oscillations, stronger than oscillations of other 
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peri ods, are believed mostly due to acoustic w aves traveling into sunspots from the surrounding ar¬ 
eas (IBeckers & Schultzll 19721: iThomasll 1985b . In the chromosp here, 3-minute oscillat ions become 
dominant and are often associated with umbral flashes (e.g., N agashima et a l. 2007). Running 
penumbral waves (RPWs) are a prominent phenomenon observed in the sunspot chromosphere, 
and they display as concentric ripples propagating from the umbra-penumbr a boundary outward 
with a typical speed of 1 0 — 25 km s" 1 (lGiovanellilll972l:IZirin & Steinlll972lt. although speeds as 
low as 4 — 12 km s -1 ( Madsen et alJl2015h and as fast as 28 — 65 km s -1 ( Kolobov & Kobanov 


2009) were also reported. The nature of the RPWs has been long debated, and in recent years 
it was suggested that these waves might represent magnet oacoustic waves that propagate upward 


from the photosphe re along inclined magnetic field lines (IBloomfield et al.ll2007l:ljess et al. 


2013 


Madsen et al. 20151. Other chromospheric observations also revealed waves propagating only in¬ 
side sunspot umbrae with a speed of 45 — 60 km s" 1 , which appeared to term inate near the umbra 
and penumbra boundaries (iKobanov & Makarchikll2004l: iLiang et al.ll201 lb . The RPWs and the 
umbral waves are thought to be independent phenomena. More recent observations from IRIS ( In¬ 
terface Region Imaging Spectrograph ) found that sunspot oscillations in the chromosphere and the 
transition region displayed the behavior of shock waves (Tian et al. 20141. 


How helioseismic waves, including both /-mode (surface gravity wave) and p-mode (acoustic 
wave), interact with sunspots has been studied to better understand the sunspots^ internal structure 
and dynamics, as well as wave properties in the presence of magnetic field. Cameron et al. (120081. 
20111 reconstructed from observations how the /-mode wave passes through a sunspot surface like 
a plane wave, and modeled numerically the sunspot’s semi-empirical internal structures through 
comparing the simulated wave with the observationally reconstructed wave. Using a similar he 


lioseismic technique to reconstruct the waves’ propagation, IZhao et alJ (201 la l and I Yang et al. 


(1201211 measured acoustic wavefunctions scattered by sunspots, and lChou et alJ (1201211 studied the 


interference fringes of acoustic waves around sunspots. From a different perspective, Zhao et al. 


(201Tb) investigated how /- and p-mode waves, originating from virtual point sources, interacted 
with a sunspot before, during, and after their encounters. Later, through numerical simulations, 
Parchevsky et al. (12011 1 studied how helioseismic waves from a point source interacted with a 
sunspot in the hope to shed light on the sunspot’s internal dynamics and magnetic field structure. 

In this Letter, em ploy ing a procedure similar to reconstruct waves traveling from virtual point 
sources (Zhao et al. 1201 lb t. we report on a detection of fast-moving waves at the photospheric 
level, propagating radially from the sunspot umbra through the penumbra to the outside. The 
observed properties of this wave are consistent with a fast magnetoacoustic wave excited beneath 
the sunspot’s surface, the wavefront of which sweeps across the photosphere. Although the cause 
of this subsurface wave source is not clear at this time, this new phenomenon may help to open a 
new window to diagnose the sunspots’ interior. We introduce the data and analysis procedure in 
Sec. 2, present the results in Sec. 3, and discuss the nature and possible cause of the wave in Sec. 4. 
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Fig. 1.— (a) Continuum intensity and (b) magnetic field of the sunspot located in AR 11312. Panel 
(a) also shows a schematic plot illustrating the data analysis scheme. Roughly, ‘O’ is the geometric 
center of the sunspot, and lines ‘QA’ and ‘05’ are along the sunspot’s radial direction. 


2. Observation and Analysis 


In this Letter, we present analysis of a sunspot located inside active region NOAA 11312 as 
one example; however, it is worth pointing out that we have obtained essentially the same results 
for several other sunspots. We believe that the results reported here are a common phenomenon 
existing in most sunspots. 


As shown in Figured! this sunspot appeared to be round, which can be presumably considered 
as axisymmetric for helioseismic purposes, and stayed relatively stable during its disk passage from 
2011 October 3 through October 17. In this study, we use the data taken by 5DO/HMI (Solar Dy - 
namics Observatory / Helioseismic and Magnetic Imager, Scherrer et al.ll2012l:ISchou et al.ll2012h 
between 00:00UT of October 8 and 23:59UT of October 12, a 5-day period during the sunspot’s 
passage from about 30° east to 30° west of the central meridian. To prepare the data for analysis, 
we track the sunspot with its local rotation rate, and remap the data into heliographic coordinates 
using Postel’s projection with the sunspot located at the center of the images. In order to assure 
that the observed phenomenon in this study is not due to artifacts or unknown calibration errors, we 
analyze four observables obtained by the SDO/HMI: Doppler velocity, continuum intensity, line- 
core intensity, and line-depth, as well as one observable from the SDO /AIA (Atmospheric Imaging 
Assembly, Lemen et al. 201 2) : 1700A line intensity. For all observables, we use running differ¬ 
ence images of the original data without applying other signal filters so that the solar convection 
gets substantially suppressed while wave signals around the sunspot are not compromised. 
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For stochastic wavefields, cross-correlating oscillation signals observed at one location with 
those observed at other locations is a powerful tool to reconstruct how waves propagate away 
from this location to oth ers. This method has been widely used in various research fields, in¬ 
cluding seismology ( e.g.. ISnieder & Wanenaan 12010: Ide Ridden l2014t) and helioseismology (e.g., 


Cameron et al 


described by 


2008; 


Zhao et al. 


Zhao et al. 

(hoilbl) . 


2011b). The analysis procedure used in this study is the same as 


but the virtual wave sources are chosen inside the sunspot. The 


analysis scheme is illustrated in Figured^: all waves originating from a virtual source ‘A’ can be re¬ 
constructed by cross-correlating the time series observed at ‘A’ with the time series observed at all 
other locations inside the field of view. Waves originating from the virtual source ‘B’ can be recon¬ 
structed similarly. However, practically, the waves constructed from a single source using a limited 
time period do not possess sufficiently high signal-to-noise ratio for further analysis. Because the 
sunspot is nearly axisymmetric, we can enhance the signal-to-noise ratio through averaging the 
wavefields obtained from sources ‘A’ and ‘B’ by rotating the wavefield from ‘B’ to overlap direc¬ 
tions of ‘OA’ and ‘OB’ and match locations ‘A’ and ‘B’. The same azimuthal averaging process 
is repeated for waves reconstructed for all other locations inside the sunspot. Moreover, when we 
average the wavefields for final results, we separate the cases for virtual sources located inside the 
umbra and penumbra in case that waves from these two different entities have different properties. 
The resultant wavefields resemble the propagation of waves originating inside the sunspot umbra 
and penumbra. However, due to that the penumbra covers a large distance range, our results from 
the sunspot penumbra smear out the differences, if any, of waves initiated from the inner, middle, 
and outer penumbra. 


3. Results 

The results for the sunspot umbra do not differ significantly from the results for the penumbra 
regarding the properties of the wave but are much noisier. Therefore, in this Letter we focus 
on presenting and discussing only the results obtained when the virtual wave sources are located 
inside the penumbra. Also, the results obtained from the different HMI and AIA observables are 
qualitatively similar, but the HMI Doppler velocity data give the best signal-to-noise ratio in our 
analysis. Figures |2] and |3] show the results obtained using the Doppler velocity data, and Figure 0 ] 
illustrates the results from the other observables. 

Figure |2] displays a few selected snapshots of the waveform obtained for the virtual source 
located inside the penumbra, and the online movie shows the entire sequence. For the time lags 
from 3.75 min to 14.25 min, a fast-moving wave is seen propagating along the sunspot’s radial 
direction from the penumbra to the outside of the sunspot. This wave is asymmetric with the 
strongest power and fastest traveling speed along the sunspot’s radial outward direction, and weak 
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Fig. 2.— Selected snapshots of the reconstructed waves propagating away from the sunspot 
penumbra, shown as foreground color images. The background black-and-white image shows 
continuum intensity of the studied region. The snapshots taken at 3.75, 5.25, 6.75, and 9.00 min 
show the fast-moving wave along the sunspot’s radial direction, and the snapshots taken at 14.25 
and 23.25 min mainly show the typical helioseismic waves expanding in all directions. An online 
movie is associated with this figure. 

or no power and slower speed for other directions. The wave signal terminates at approximately 
30 Mm away from the sunspot center, or about 15 Mm beyond the sunspot boundary. This fast- 
moving wave is clearly distinguishable from the usual helioseismic waves, which can be seen in 
the snapshots at 14.25 min and 23.25 min. Actually, in our analysis, the helioseismic wave is first 
identified only at 12.00 min. 

To better quantify the properties of the fast-moving wave, we obtain the time-distance relation 
of the wave by stacking together the cross-correlation functions for different travel-time lags. Since 
the wave travels preferentially along the radially outward direction, only the time-distance diagram 
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Fig. 3.— (a) Time-distance diagram obtained along the sunspot’s radially outward direction, (b) 
Same as (a), but with two white boxes delimiting the areas of the fast-moving wave (lower box) 
and helioseismic waves (upper box). The green line shows a linear fitting of the fast-moving wave, 
corresponding to a speed of 45.3 km s -1 . (c) Time-distance diagram from a quiet-Sun region is 
shown for comparison, (d) Power-spectrum diagram for the fast-moving wave. The green line 
corresponds to the same phase speed as the line in panel (b). 


obtained along this direction is shown in Figure|3j There are two wave branches visible in the time- 
distance diagram (Figure [3^-b), with the lower branch corresponding to the fast-moving wave, and 
the upper branch corresponding to the usual helioseismic waves. By comparing this time-distance 
diagram with the one obtained from a quiet-Sun region following the same analysis procedure 
(Figure [3]:), we find that the helioseismic waves from the penumbral sources are largely suppressed 
for short travel distances, and that the fast-moving wave is a phenomenon only associated with the 
sunspot. 

The fast-moving wave appears to propagate with a nearly constant speed and shows little 
dispersion, i.e., its group velocity is similar to its phase velocity. We use a linear fitting to estimate 
its apparent phase velocity and get a speed of 45.3 ±1.7 km s -1 , substantially faster than the speed 
of fast magnetoacoustic wave or Alfven wave, which is an order of 10 km s -1 in the photosphere 
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of sunspots. Taking the time-distance diagram of the fast-moving wave (the lower white box in 
Figure [3 ]d), we calculate the power spectrum of this wave (Figure [3}i). Due to the limited spatial 
and temporal scale of the wave, aliasing is hard to avoid when performing Fourier transform and 
calculating the power spectrum. The k — v power-spectrum diagram shows a dominant power in 
the range of 2.5 — 4.0 mHz, indicating that the oscillation frequency of this fast-moving wave falls 
into the category of 5-minute oscillation, same as the typical acoustic waves in the photosphere. 
The dominant power shows a phase velocity close to 45 km s _ 1 , in agreement with the fitting result 
in the space-time domain. 

The fast-moving wave is also detected in the analysis results from the other HMI and AIA 
observables (Figure 3k-d), although the signal-to-noise ratio varies among these observables. The 
properties of the observed fast-moving wave are also slightly different in the characteristic fre¬ 
quency and phase velocity, probably because different observables correspond to different atmo¬ 
spheric heights and the wave has different properties at different heights. 

As a matter of fact, although the fast-moving wave is not apparent in the HMI Doppler 
and continuum intensity data without applying the cross-correlation calculations, part of the fast- 
moving wave is directly visible in the HMI line-core and line-depth data. Figure [4je displays a 
time-space plot of the line-core data with a horizontal cut through the sunspot center, in which 
one can see herringbone patterns in the inner penumbra, i.e., about 7-12 Mm away from the um¬ 
bra center delimited in the dashed lines in Figure 0k . The herringbone patterns have a propagation 
speed of approximately 45 km s -1 , consistent with the estimate obtained from the cross-correlation 
calculations. Between 12-20 Mm, the outer penumbra is dominated by slow-moving waves with 
an apparent speed of 6 km s -1 . These waves seem to be an extension of some of the fast-moving 
waves, but they do not appear in the cross-correlation diagram (Figure HJd), implying these waves 
do not share common sources with the fast-moving waves. Although not clearly visible in Figured 
beyond 12 Mm, the fast-moving waves can still be picked up by the cross-correlation analysis up 
to 30 Mm from the sunspot center. What is the relation between the fast-moving wave in the inner 
penumbra and the slow-moving wave in the outer penumbra is worth further studying. 


4. Discussion 

Through cross-correlating stochastic oscillation signals observed at various locations inside 
and near a sunspot, we reconstruct waves propagating away from virtual point sources located 
inside the sunspot. It is found that after a virtual source is initiated and before typical helioseis- 
mic waves are visible, a fast-moving wave appears propagating outward along the sunspot’s radial 
direction, with an apparent phase velocity of about 45 km s _1 and a characteristic frequency be¬ 
tween 2.5 and 4.0 mHz. This wave starts from the sunspot umbra and terminates at about 15 Mm 
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Fig. 4.— Time-distance diagrams obtained in the same way as for the diagram shown in Figured, 
but using data of (a) HMI continuum intensity, (b) HMI line-core intensity, (c) HMI line depth, and 
(d) AIA 1700A. (e) A time - space plot obtained with a horizontal cut through the sunspot center, 
displayed in the HMI line-core intensity data. The inclined dashed lines indicate a propagation 
speed of 45.3 km s -1 . The waves of interest are confined between the white dashed lines, about 
7 — 12 Mm from the sunspot center. Other wave patterns are visible beyond this range but are not 
related to the phenomenon discussed in the text. 

beyond the sunspot boundary. In the inner penumbra, for about 1/4 of the whole traveling range, 
the wave is visible in the inner penumbra in the HMI line-core intensity data without applying 
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Fig. 5.— (a) Schematic plot showing a scenario of that a subsurface disturbance generates a 
helioseismic wave sweeping across the photosphere, forming the fast-moving wave observed in 
the photosphere, (b) Vertical snapshot showing wave propagating away from the disturbance, 
located at 2 Mm away from the sunspot’s central axis at the depth of 5 Mm. The snapshot is 
taken at 7.5 min after the wave is excited. The wave shown in this figure is obtained from solving 
the ray-theory MHD equations using a realistic sunspot model for the frequency of 3 mHz. (c) 
Time-distance diagram seen in the photosphere for this wave. For comparison, a white dashed 
line, representing a speed of 45.3 km s -1 is plotted. 

cross-correlation calculations. For the rest of the wave traveling range, and for HMI Doppler and 
continuum intensity data, the wave is only visible after cross-correlation analysis is applied. 
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What is the nature of this fast-moving wave, what causes it, and how to explain its fast speed 
are interesting questions. Although Figure [2] (and the online movie) seems to show the wave is a 
traveling wave on the surface, it is physically implausible to explain its propagation speed, which is 
substantially faster than the local Alfven wave and fast magnetoacoustic wave. The characteristic 
frequency of the wave implies that this wave falls into the category of acoustic waves, and the time- 
distance relation of the wave is a key factor to determine the nature of the wave. To explain this 
wave, we conjecture a disturbance occurring at approximately 5 Mm beneath the sunspot surface 
and 2 Mm away from the sunspot’s central axis (see Figure |5k). Acoustic waves are excited by 
this disturbance and expand toward all directions. Due to the stratified structure of the sound- 
speed profile as a function of depth, and the existence of the magnetic field that will alter, albeit 
slightly, the sound-speed profile both horizontally and vertically, the wavefront does not expand 
symmetrically toward all directions. The wavefront touches different locations of the photosphere 
at different times, forming effectively a fast-expanding ellipse at the surface due to the asymmetric 
propagation speed. Numerous disturbance sources are located below the sunspot in a wide range 
of depths and distances to the central axis, and cross-correlation analysis eventually gives us a 
wave as seen in Figure |2l That is, it is likely that what we see at the surface are wavefronts 
of the waves, which are excited below sunspots, sweeping across the photosphere. Employing 
a realistic magnetohydrodynamic sunspot model (Remoel et al. 2009) and following the ray-path 
approximation with magnetic field (Moradi & Cally 2008), we calculate how a magnetoacoustic 
wave, at a frequency of 3.0 mHz, expand from the deep source and how its wavefront sweeps 
across the photosphere (see Figure O). The time-distance relation from this calculation shows an 
approximately 40 km s -1 surface speed of the wavefront when it is 10 Mm horizontally away from 
the source on the surface. The difference in the apparent speeds between the model and observation 
may be due to the relatively weak magnetic field used in the model, or due to different depths of 
disturbance sources. We only present here our calculations of 3.0 mHz, but we also recognize 
that the subsurface disturbances can excite waves with very broad frequency band. The frequency- 
dependent power distribution at the subsurface locations is determined by the local atmospheric 
properties, and is expected to be different from the power distribution observed at the photosphere. 

However, despite that this conjecture is promising in explaining the observations, there are 
clear difficulties. One difficulty is that it is unclear what causes those disturbances beneath sunspots? 
Some past helioseismic analyses indicated that the depth of 5 Mm was approximately where the 
sunspot’s downdraft meets the upward flow (e.g., Zhao et al. 2010), and plasma collisions near 
this depth may cause disturbances that generate the magnetoacoustic waves. Other than that, cur¬ 
rently there is no other evidence as we know indicating that this area may have plentiful wave 
sources. The other difficulty is, the magnetohydrodynamic processes between sunspots’ subsur¬ 
face and surface are rather complicated, and our equations prescribing the ray paths do not take 
into consideration of the subsurface flow field, the alteration of wave eignfunctions in the presence 
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of strong magnetic field, and maybe other factors. The inclined magnetic field and fast flows in 
the penumbra further complicate the estimate of travel speed of the waves. Therefore, we do not 
expect that the picture given in Figure |5] fully represents all details of the wave propagation through 
the sunspot region, but it is a useful first attempt to explain the observed phenomenon. Meanwhile, 
we cannot completely rule out other possible explanations, e.g., what we observe may be waves 
reflected back from the upper chromosphere (Cally & Moradi 20131. 


It is likely that the chromospheric RPWs (see §1) and the photospheric fast-moving waves 
have a common cause. When a disturbance occurs about 5 Mm below the sunspot surface near 
the sunspot’s central axis, a magnetoacoustic wave sweeps across the photosphere in the sunspot 
and its immediate vicinity, and a fast-moving wavefront is detected. When the wave propagates 
along the inclined magnetic field from the subsurface to the photosphere and above, the RPW 
is observed in the chromosphere as the transverse projection of the magnetoacoustic wave (e.g., 
Bloomfield et al. 2007 1. 


This newly detected fast-moving wave, although may not be a new type of wave, likely carries 
rich information from sunspots’ subsurface area, and may help to open a new window to diagnose 
sunspots’ internal structure and dynamics. 


SDO is a NASA mission, and HMI project is supported by NASA contract NAS5-02139 to 
Stanford University. 
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